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simulation was performed using SU(2) lattice gauge theory in the maximal Abelian 
gauge. We find that the monopole current and electric flux distributions are quite 
different below and above the finite temperature deconfining phase transition point 
and agree with predictions of the Ginzburg-Landau theory. 
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I. INTRODUCTION 



It is widely accepted that quark confinement is due to color-electric flux tube formation. 
In this scenario the QCD vacuum behaves like a dual superconductor. That is, magnetic 
monopoles in the QCD vacuum respond to a color-electric field by generating screening 
currents to confine the electric flux into a narrow tube dual to the Abrikosov vortex arising 
from the Meissner effect. Some direct evidence for this scenario has been reported for U(l) 
LGT and SU(N) lattice gauge theories (LGT) simulations ]]]-§]. It is found that monopoles 
in U(l) and the [/(l)^ -1 abelian projection of SU(N) LGT in the maximal Abelian gauge @ 
react to the electric flux from a static qq pair by producing a solenoidal current distribution 
which confines the flux to a narrow tube. 

It is clear that the solenoidal magnetic monopole currents are associated with the con- 
fined flux tube. At finite temperatures the flux tube is expected to exist only in the low- 
temperature confined phase, but not in the high-temperature unconfined phase. Previous 
studies addressed the issues of monopole current distributions and the structures of 

the electric flux tube at zero temperature. Questions remain as to the monopole current 
distributions at finite temperature, especially in the unconfined phase. In this paper we 
shall present our observations of monopole current and electric flux distributions at finite 
temperature. This study is a direct extension of the earlier work in SU(2) LGT. 



II. SIMULATIONS 

Our simulations were performed on a lattice of the size 4 x 17 2 x 19 with skew-periodic 
boundary conditions. Each link from point s in the /i direction carries an £77(2) element 
U (s, /i), and the plaquette operator Up(s) is formed in the usual fashion as a directed product 
of link variables, where the label P represents the orientation of the plaquette, (/i, v). The 
standard Wilson action is used in our calculations 

S(U)=(3J2(1-I^U P ). (1) 

S,P 
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We generated lattice configurations with the distribution exp(— S) using a combination of 
the standard Metropolis algorithm || and overrelaxation ||. Simulations were performed 
for (3 = 2.25, 2.28 and 2.40. In each case the initial 3000 sweeps were used to thermalize the 
system, then one measurement was made every 50 sweeps. 

Each measurement was performed on a lattice configuration in the maximal Abelian 
gauge H, which can be achived by finding the gauge transformation that maximized the 
quantity |10] 



R = Y,TA°zU(s^)a 3 U^s^)}. (2) 

This is equivalent to diagonalizing 

= £[t/(^y 3 C/W) + U"^(s — n, n)a 3 U(s — /x, //)]. (3) 

at each site s. To measure the efficacy of gauge fixing in the simulations we used the lattice 
sum of the magnitude of the off-diagonal component of X(s), \Z\ 2 = J2 S |A(s)i 2 | 2 . Typically 
we needed about 620 gauge fixing sweeps to attain \Z\ 2 ps 10~ 5 /site for each (3 value. Three 
gauge fixing methods were used: (1) generating and accepting random local changes only if 
R increased, (2) locally maximizing R exactly at alternate sites, (3) applying overrelaxation 
using the square of the gauge transformation of method (2) to sample configurations better. 



III. MEASUREMENTS 

After gauge fixing, the Abelian U(l) link variable is given |H| by the phase of the diagonal 
component of the SU(2) link variable, u(s,fi) = U(s, fi) u /\U(s, /x)n|. We construct the 
Abelian plaquette variable p^ u {.s) from m(s,/i) in the usual fashion. In this study we use 
Abelian Polyakov loops to replace Abelian Wilson loop used in ref. 0,0], in order to study 
the finite temperature case. The Abelian Polyakov loop closed in the time direction can be 
contructed as 

Nt 

p ab (s) = n«(v), (4) 

r=l 

3 



where N t is the temporal size of the lattice, N t x iV s x N z , specifically N t = 4 for the lattice 
4 x 17 2 x 19. 



The magnetic monopoles are identified using the DeGrand-Toussaint construction ||TT 
For each Abelian plaquette p^s) we decompose its phase angle 9 P = arg{p IJ , v {s)) (— 4n < 
9 P < 4tt) into two parts, 

P = 9 P + 2ixn pi (5) 

with 7i < Op < 7r, where n p describes the number of Dirac strings and V the electromagnetic 
flux through the plaquette RTl |. The magnetic monopole current density can be constructed 



as 12 



K^s) = t^ pu [O p {s + i/, pa) - p (s,pa)]/4ir. (6) 

It is convenient to associate the monopole current density K^(s) in each three- volume with 
a link on the dual lattice 0, making world lines that define a conserved current density Jm- 
To isolate the solenoidal monopole currents we construct the operator for the line integral 
of Jm around a dual plaquette, V X J ]y[ , clS described in ref. H. 

In our measurements we chose two Abelian Polyakov loops correlated along the N z 
direction, (P a b(0)P^ b (R)), with the separation R = 3a, 4a, 5a and 6a, which represent the 
static qq pair. The electric field operator is given by a 2 Si = Impj 4 , where a is the lattice 
spacing M. Then the electric field in the N z direction can be calculated as, 



(P afe (0)P a t b (it!)Im(^(x ))) 
a*e(P ab (0)Pl(R)) 



with e the Abelian electric charge and p z t the Abelian plaquette in the z — t plane. The 
solenoidal monopole current distribution can be calculated from the correlation [H , 



2n(P a b (0)Pj b (R)(V x J M (x)) z ) 
a*e(P ab {Q)Pi b {R)) 



(v x Jm {r^)) z = ^ l ^ t 7;;^ t ;: MW ^ ( 8) 



where the solenoidal monopole current operator, (V x Jm)z, is constructed as the dual 
plaquette in the x — y plane. In practice, the expressions for (E Z (R, x)) and (V x Jm(R, x ))z 



can be simplified further. If we denote 6ppt as the phase angle of the product of two Abelian 
Polyakov loops, (P o &(0)Pj 6 (i£)), Eq. (|7|) can be written as, 

( g ,(j;,x))= <t"y , o) 

a z e{cosUp P i) 

where 9 p is the phase angle of the Abelian plaquette p zt . Also, Eq. (||) becomes 

(vxj M (fi,x))^ 2 ^ s '"V; (v ; J r (x)>j> - < io > 

a 4 e(cosfc'ppt) 



IV. NUMERICAL RESULTS 

As we choose the lattice extent N t of the 4-dimensional lattice, N t x N% x N z , as the time 
direction, we can simulate the finite temperature case on the lattice with the temperature 
T = l/N t a(P), where the lattice spacing a is a function of (3. For SU(2) lattice gauge theory 
the finite temperature deconfining phase transition point has been studied extensively. One 



recent simulation result is, for N t = 4 the transition point f3 c is given as, 

P c = 2.2986 ± 0.0006. (11) 

Below (3 C ((3 < f3 c ) the system is in the confined phase, but above f3 c it is in the unconfined 
phase. 

We measured the electric flux (E z (R,x)) and the solenoidal monopole current (V x 
Jm(R, x )) 2 on the transverse plane midway between the qq pair for several (3 values, (3 = 2.25, 
2.28 and 2.40. Two (3 values (2.25 and 2.28) are below the transition point (3 C given by 
Eq. (|i~l~D , one (3 value (2.40) is above the (3 C . In the following we proceed to discuss our results 
in the two cases, which correspond to the confined and the unconfined phases separately. 

A. Results in the confined phase 

We measured the electric flux and the solenoidal monopole current distributions of a 
qq pair in the confined phase: (3 = 2.25 and 2.28, In each case we accumulated about 800 
measurements. 



Fig. p] shows the electric field E z distribution for (3 = 2.28, with the qq separation R = 3a, 
4a, 5a and 6a separately. This displays the strength of E z versus the distance r from the 
qq axis on the transverse plane midway between the qq pair. For small quark separations, 
e.g. R = 3a and 4a, as shown in Fig. |l|(a) and (b), our results are consistent with those of 
the earlier works, which were obtained from simulations using Abelian Wilson loops of sizes 
3 x 3 in ref. || and 3x3 and 5 x 5 in ref. f|. In this study we observe clear signals for E z 
at R = 5a, 6a as shown in Fig. |l|(c) and (d). From Fig. [I] one can see that the electric flux 
E z has the peak value on the qq axis (r = 0), then it decreases with the off-axis distance 
r rapidly. As the qq separation R increases the E z values also decrease as we expected. 
Our data also shows that the electric flux E z basically vanishes within errors at the off-axis 
distance r pa 5a for all R values in Fig. |l|(a)-(d). This agrees with the expectation that 
electric flux is squeezed into a narrow tube in the confined phase. 

In Fig. ^|we show the monopole current — (V x Jm)z as a function of the off-axis distance 
r for f3 = 2.28, with R = 3a and 4a. Again our data shows that the curl of the monopol 
current — (V x Jm)z has the same functional form as those of the earlier works for SU(2) 
LGT Fig. only displays the results at small quark separations, e.g. R = 3a and 4a 

since for large quark separations, e.g. R = 5a and 6a the noise is too large. The large signals 
shown in Fig. ^| confirm that the solenoidal monopole current surrounding the electric flux 
tube is large in the finite temperature confined phase. For (3 = 2.25 the results are similar 
to the case of (5 = 2.28. 

To clarify the significance of these profiles, we wish to summarize briefly the results of ref. 
Jl|fJ]. Picture for the moment the points in Fig. [| at r/a = 1, 1.4 to be negative and further 
that the curve matchs, up to a factor, the negative of the electric field profile for r/a ^ 0. 
This is precisely what we found for U(l), ref. 0, and is the prediction of the London theory. 
The linear combination of E z — A 2 (V x Jm)z is defined as the fluxoid, where A is the London 
penetration depth. For this geometry the fluxoid gives a lattice delta function at the origin, 
with the coefficient equal to the total electric flux in the vortex. 

The departure of Fig. from this picture is fully explained by the more general Ginzburg- 

6 



Landau theory in which the order parameter for superconductivity turns on at a boundary 
over a finite distance (called the coherence length) rather than abruptly at the boundary as 
in the above case. The coherence length in the Fig. ^| is the distance at which the profile 
deviates from the London theory which is r/a =~ 2.0. The London penetration depth is 
independently determined by the exponential fall-off of E z . 

B. Results in the unconfined phase 

We also studied the electric flux and monopole current distributions in the unconfined 
phase. The data were measured at ft — 2.40, and 410 measurements were accumulated. 

Fig. |H shows the electric flux E z as a function of the off-axis distance r for (3 = 2.40, 
with the qq separation R = 3a, 4a, 5a and 6a separately. The significant difference between 
Fig. [3] and Fig. |1] is that the E z flux values for f3 = 2.40 approach zero very slowly, even 
at large off-axis distances (e.g. r = 6a) E z flux values still do not vanish within errors, as 
shown clearly in Fig. |3](c) and (d). However, Fig. [I] shows that the E z flux data for f3 = 2.28 
decrease very quickly and vanish at r « 5a within errors. This implies that the electric flux 
spreads out from the qq axis in the unconfined phase {(3 = 2.40), it is more Coulomb-like 
than the flux distribution in the confined phase (e.g. (3 = 2.28). 

Fig. f| shows the corresponding monopole current distributions —(V x Jm)z for /3 = 2.40, 
with R = 3a and 4a. From this figure one can see that almost all data vanish within errors 
except one small value on the qq axis (r = 0) for each case. Comparing with the monopole 
current distributions in the confined phase (/3 = 2.28), as shown in Fig. § it is clear that 
the solenoidal monopol current — (V x Jm) z becomes very small in the unconfined phase, as 
expected. 
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V. CONCLUSIONS 



We extend the earlier studies |2||| on solenoidal monopole current and electric flux distri- 
butions to the finite temperature case. We find that both the electric flux and the monopole 
current have quite different behaviors in the finite temperature confined and unconfined 
phases. In the confined phase our data are consistent with the results of earlier works , 
which correspond to the zero-temperature case. We find that in the confined phase the 
electric flux is confined in a narrow tube, and the solenoidal monopole currents are large. 
However, in the unconfined phase our data show clearly that the electric flux is not con- 
fined, and the solenoidal monopole current is vanishingly small. This presents the significant 
difference between the two phases. This study confirms the point that the electric flux is 
repelled by the surrounding solenoidal monopole currents, and is squeezed into a narrow 
tube. In the unconfined phase solenoidal monopole current vanishes, then the electric flux 
is unconfined. 
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FIGURES 

FIG. 1. Profile of the electric field E z as a function of the transverse distance r from the qq 
axis, with the qq separation (a). R = 3a, (b) R = 4a, (c) R = 5a and (d) R = 6a. The data were 
measured on the lattice 4 x 17 2 x 19 with (3 = 2.28. The values of E z are measured in the lattice 
unit l/a 2 e. The dashed lines are just to guide the eye. 

FIG. 2. Profile of the solenoidal monopol current — (V x Jm)z as a function of r, with the qq 
separation (a). R = 3a and (b) R = 4a. The data were measured on the lattice with (3 = 2.28. 
The values of —(V x Jm)z are measured in the lattice unit l/a 4 e. 

FIG. 3. Profile of E z vs. r in the unconfined phase with (3 = 2 AO. The qq separation is (a). 
R = 3a, (b) R = 4a, (c) R = 5a and (d) R = 6a. The values of E z are measured in the lattice unit 
l/a 2 e. 

FIG. 4. Profile of (V x Jm)z vs. r in the unconfined phase with (3 = 2.40. The qq separation 
is (a). R = 3a and (b) R = 4a. The values of — (V x Jm)z are measured in the lattice unit l/a 4 e. 
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